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We numerically investigate the three-dimensional 0(6) model on 12 3 to 120 3 lat- 
tices. From Binder's cumulant at vanishing magnetic field we obtain the critical 
coupling J c = 1.42865(5) and verify this value with the x 2 -method. The universal 
value of Binder's cumulant at this point is g r (Jc) = —1.94456(10). At the critical 
coupling we find the critical exponents v = 0.818(5), /3 = 0.425(2) and 7 = 1.604(6) 
from a finite size scaling analysis. We also determine the finite-size-scaling function 
on the critical line and the equation of state. Our 0(6)-result for the equation of 
state is compared to the Ising, 0(2) and 0(4) results. 



1. Introduction 



Our goal is the determination of universal quantities of the three- 
d ' dimensional 0(6)-invariant nonlinear er-model, especially critical exponents 

and scaling functions. These quantities shall eventually be compared to 
those of 2-flavour staggered QCD with adjoint fermions. The relevant sym- 
metry of this model is SU(4), so it should be in the same universality class 
as the 0(6) spin model 1 . The Hamiltonian of our model is defined as 

PH = -J £ $ x -$v - H- (!) 

<x,y> x 

Here x and y are the nearest-neighbour sites on a three-dimensional hy- 
percubic lattice, (j> x is a 6-component unit vector at site x and H is the 
external magnetic field. The coupling constant J is considered as inverse 
temperature, therefore J = 1/T. At H = \H\ ^ the order parameter, the 
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magnetisation M, is defined by 

M = (d>W), (2) 

where is the longitudinal (parallel to H) component of 4>. At H = this 
quantity vanishes for all couplings on the lattice, but we can use 

M = (\$\) (3) 

as an approximate order parameter. We also measure the susceptibility x 
and Binder's cumulant g r , defined by 

X = V({ 4? ) — M 2 ), (4) 
* = ^-3. (5) 

(02)2 

We have simulated the model with a Wolff-Cluster- Algorithm, using a 
ghostspin to emulate the external magnetic field. We used hypercubic lat- 
tices with periodic boundaries and linear extensions L between 12 and 120. 
Near the critical coupling J c we have made up to 200000 measurements per 
data point, elsewhere about 20000 measurements. Between the measure- 
ments we have used up to 400 cluster updates in the broken phase and up 
to 1500 updates in the symmetric phase. 



2. Determination of the critical coupling 

At H = the Binder cumulant can be described by the finite-size-scaling 
function 

g r = Q g {tL 1 ^;L- u ), (6) 

where t = T 7 r T " is the reduced temperature and u> is the leading irrelevant 
exponent. Therefore g r should be independent of L at the critical point 
t = apart from corrections due to irrelevant scaling fields. The left plot 
in Fig. 1 shows our reweighted data, with the dotted lines representing the 
jackknifc error corridors. Although he intersection points coincide within 
the errors, there are some small corrections. As the value of u is unknown, 
we used Binder's approximation 2 

— = 1 + — (7) 

J ip J c log b 

to extrapolate the intersection points Ji P {L,b) of two lattices with sizes L 
and L' = bL. Fitting the results for the lattice sizes L = 12, 16, 20, 24, 30, 36 
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Figure 1. The left plot shows the binder cumulant for the lattices with size L = 
12,16,20,24,30,36,48,60,72. The large lattices have a steeper slope then the small lat- 
tices. The right plot shows the finite-size-scaling function at the critical coupling. 



0.699960(14) 



1.42865(5). 



(8) 



to a constant value we find 
1 

T c 

We have checked this result with the x 2 -method described in Ref. 3. From 
this value we derive the universal quantity g r (J c ) = —1.94456(10). 
The result for J c is comparable with the result J c = 1.42895(6) of Butera 
and Comi 4 using a high temperature expansion, but their errors seem to 
be underestimated. 



3. The critical exponents /3, 7 and v 

For the determination of the critical exponents we use the following scaling 
relations at t = H = 0: 

M - L-V" (ao + a 3 L~ u ) (9) 
X - Vl" (6 + b 3 L-») (10) 

^ = L x l v (do + d 3 L- w ) (11) 

Fitting these relations to our data at lattices in the range L = 12 — 72 
we get P/v = 0.519(2) and 7/1/ = 1.961(3) from M and X - The error of 
these quantities includes a variation of 10 between 0.5 and 1.0. From the 
derivative of g r we get 1/v = 1.223(5). Here we neglect the correction term, 
because it is zero within the errors. From these results we can calculate the 
critical exponents (3 = 0.425(2), 7 = 1.604(6) and v = 0.818(5) . They fulfil 
the hyperscaling relations and are in complete agreement with the results 
of Butera and Comi. 4 



February 1, 2008 4:3 WSPC/Trim Size: 9in x 6in for Proceedings 



proc'heidelberg 



4 

4. The equation of state 

In the vicinity of T c the critical behaviour can be described by the universal 
equation of state. It can be written in the form 

M = h 1/5 f G {z'), z' = t'/h 1/fi5 (12) 

t' and h are the normalized reduced temperature t' = (T — T c )/To and 
magnetic field h = H/H n with the normalization conditions /g(0) = 1 and 
/g(V) = (— z')P as z' — > — oo. To get the normalization constants Ho and 
To we have to determine the critical amplitudes of the magnetisation on the 
critical isotherm and on the coexistence line. The infinite volume behaviour 
of M is given by 

M(T C ,H) = doH 1 ' 5 ^ + d\H" v *) (13) 

including the leading correction term. From hyperscaling relations we get 
for the exponents 5 = 4.780(22) and v c = 0.4031(24). With this ansatz, 
using only the data of the largest lattice at each value of H, we get for the 
critical amplitude d c = 0.642(1), which translates to H = d~ s = 8.3(1). 
In the broken phase the magnetisation is described by the ansatz 

M(T < T c , H) = M(T, 0) + Cl (T)H^ 2 + c 2 (T)H. (14) 

The iJ 1 / 2 term is due to the Goldstone effect. As we have already seen 
with the 0(2) 5 and 0(4) 6 spin model, this ansatz fits very well and allows 
the extrapolation of data at nonzero H to H = 0. Fitting the results for 
M(T, 0) for several couplings (J = 1.45,1.47, 1.5,1.55 and 1.6) to the ansatz 

M(T < T c , 0) = B (T c - Tf [1 + h (T c - Ty v + b 2 (T c - T)} (15) 

we get B = 1.22(1), which leads to T = B^ 1 ^ 3 = 0.63(1). 
The left plot in Fig. 2 shows the scaling function /g(^')- I n the broken 
phase (z' < 0) the solid lines are reweighted data from different J-values. 
There are visible corrections, so we use the ansatz 

Mh- 1 ' 6 = f G (z>) + h^f^(z') + h 2 ^f^(z') (16) 

to obtain the universal value of /g(z') (dashed line). In the symmetric 
phase the corrections are negligible. In the right plot of Fig. 2 we show a 
comparison of the 0(6) scaling functions to those of some other O(N) mod- 
els, which have been determined in Refs. 7 (Ising), 5 (0(2)) and 6 (0(4)). 
As one can see, the functions get steeper with increasing N. These func- 
tions can eventually be used to determine the universality class of staggered 
2-flavour QCD both with fundamental or with adjoint fermions. But until 
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Figure 2. The equation of state of 0(6) (left plot) and the comparison of the equations 
of state of O(l) (Ising), 0(2), 0(4) and 0(6) (right plot). 

now QCD data has too low statistics to distinguish between these scaling 
functions, and also the used lattice sizes are too small in QCD to be sure 
that one has reached the infinite volume limit. 

5. Finite-size-scaling function at J c 

Because of the small lattice sizes used in QCD, finite-size-scaling functions 
are a better tool for the comparison of spin models to QCD. At T c and for 
small H the different lattices should scale like 

M(T C ,H,L) = L~®l v Q M (HL f3S ^) (17) 

where Qm is a universal scaling function. Our results are shown in the 
right plot in Fig. 1. The solid line shows the asymptotic behaviour in the 
limit L — > oo 

Qm(z) = QM = d c z x l\ (18) 
which is observable for z = HL 0S /' / >AO. 
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